Critical Tunneling Currents in Quantum Hall Superfluids: 
Pseudospin- Transfer Torque Theory 
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At total filling factor v = 1 quantum Hall bilayers can have an ordered ground state with spon- 
taneous interlayer phase coherence. The ordered state is signaled experimentally by dramatically 
enhanced interlayer tunnel conductances at low bias voltages; at larger bias voltages inter-layer 
currents are similar to those of the disordered state. We associate this change in behavior with 
the existence of a critical current beyond which static inter-layer phase differences cannot be main- 
tained, and examine the dependence of this critical current on sample geometry, phase stiffness, and 
the coherent tunneling energy density. Our analysis is based in part on analogies between coherent 
bilayer behavior and spin-transfer torque physics in metallic ferromagnets. Comparison with recent 
experiments suggests that disorder can dramatically suppress critical currents. 
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I. INTRODUCTION 

At Landau level filling factor v = 1 bilayer two- 
dimensional electron systems in the quantum Hall regime 
can have broken symmetry ground states^"— with sponta- 
neous inter-layer phase coherence. These ordered states 
can be viewed as excitonic superfluids,— or as XY 
pseudospin ferromagnets^ in which the pseudospin is 
formed from the two-valued which layer quantum de- 
gree of freedom. The most robust experimental signa- 
ture of these states, a vastly enhanced inter-layer tunnel 
conductance^""-^ at small bias voltages, is still poorly un- 
derstood from a quantitative point of view. 

Two types of ideas, which differ most essentially in how 
the bias voltage is introduced in the theory, have been ex- 
plored in an effort to understand the height and width of 
the tunnel conductance peak. In one approacbi^rJ^ the 
bias voltage is introduced as an effective magnetic field, 
uniform across the bilayer, which induces pseudospin pre- 
cession around the z axis, driving the interlayer phase 
difference at a steady rate and inducing a purely oscil- 
lating inter- layer current. When the microscopic inter- 
layer tunneling amplitude is treated as a perturbation, 
thermal and disorder fluctuations of the condensate are 
then responsible for a finite dc conductance peak. We 
refer to this type of theory below as the weak-coupling 
theory of the tunneling anomaly. Weak-coupling theo- 
ries prediclji^""— splitting of the zero voltage tunnel con- 
ductance peak into separate finite voltage peaks in the 
presence of a magnetic field component parallel to the 
two-dimensional layers, an effect that is not^ seen ex- 
perimentally. The second type of transport theorjJ^""— 
is formulated in terms of local chemical potentials of 
fermionic quasiparticles which may be altered by the or- 
dered state condensate but are still responsible for charge 
conduction. In this type of theory the tunnel conduc- 
tance is finite even in the absence of disorder and ther- 
mal fluctuations because charge has to be driven between 
normal-metal source and drain contacts. The resistance 



generally dependsiSii^l on how the fermionic degrees of 
freedom which carry charge between leads are influenced 
by order parameter and disorder configurations and can- 
not be described in terms of condensate dynamics alone. 
In the second approach the width in voltage of the con- 
ductance peak is simply equal to the product of this re- 
sistance and the maximum current between source and 
drain at which the order parameter can maintain a time- 
independent steady-state value. 

In this article we expand on the second type of the- 
ory of the conductance peak, using pseudospin-transfer 
torque ideas^i^^ borrowed from recent ferromagnetic- 
metal spintronics literature to model the influence of the 
transport current on the pseudospin magnetization. We 
find that the critical current depends in general on de- 
tails of the sample geometry and on how disorder and lo- 
calization physics influence transport inside the system. 
Generally speaking however, the critical current is pro- 
portional to system area and to the inter-layer tunnel- 
ing amplitude when the condensate's Josephson length is 
larger than the system perimeter, and proportional to the 
system perimeter and to the square root of the inter-layer 
tunneling amplitude when it is shorter. 

Our paper is organized as follows. In Section II we in- 
troduce the pseudospin-transfer torque theory of order- 
parameter dynamics in a bilayer quantum Hall ferromag- 
net. We use this theory in Section HI to discuss critical 
current values from a qualitative point of view. In Section 
IV we report on a series of numerical studies which take 
into account the two-dimensional nature of the systems 
of interest and the edge dominated current-paths typical 
of strong magnetic fields. The pseudospin transfer torque 
theory enables us to assess the influence of sample geom- 
etry on critical currents. Finally in Section V we discuss 
the significance of our findings in relation to recent exper- 
iments. We find that experimental critical currents are 
several orders of magnitude smaller than theoretical ones 
and argue that vortex-like disorder-induced pseudospin 
textures must be largely responsible for this discrepancy. 
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We propose new experimental studies which can test our 
ideas and thereby achieve progress toward a quantitative 
theory of the spontaneous coherence tunnehng anomaly. 



II. PSEUDOSPIN MAGNETISM AND THE 
LANDAU-LIFSHITZ-SLONCZEWSKI EQUATION 

In the lattice model^i of = 1 bilayer systems, the 
local which layer degree of freedom can be expressed as 
a pseudospin defined by the operator 



(1) 



where i is the site index, a is the layer index and f is the 
Pauli matrix vector. The pseudospin Hamiltonian has 
the forro^i 



where Hij — {i,a; j,a'\Vcoi\i,o'; ja') is the direct 
Coulomb interaction associated with the z pseudospin 
component {i.e. with charge transfer between layers), 
Ffj = {i,a;j,a\Vcoi\i,o';ja), is the exchange interac- 
tion between orbitals in the same layer and Ffj — 
(«, (t; j, ct|T4oi|*, cr; ja) is the exchange interaction between 
orbitals located in different layers. Since Hij is gener- 
ally larger than Ffj , the classical ground state is an easy- 
plane pseudospin ferromagnet with a hard z axis. In the 
limit of smooth textures the pseudospin energy functional 
has the formi 
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where m = {mx,my,mz} is the local pseudospin direc- 
tion. The parameters which appear in this expression are 
the anisotropy parameter /3 > 0, the pseudospin stiffness 
(or equivalently the exciton superfluid density) ps, the 
splitting between symmetric and antisymmetric single- 
particle bilayer states due to interlayer tunneling At, and 
the pseudospin density n. 

The mean- field-theory pseudospin ferromagnet 1 con- 
sists of a full Landau level of electrons in identical phase 
coherent bilayer states. It follows that the mean- field- 
theory pseudospin density n is equal to the full Landau 
level density, {2TrP)-^. (Here I = {hc/eBf/'^, where B 
is the magnetic field strength, is the magnetic length.) 
Mean-field-theory can also used^ to find explicit expres- 
sions for Ps and /3. In practice the values of these three 
parameters are modified^^ by quantum and thermal fluc- 
tuations. The fourth parameter Aj is exponentially sen- 
sitive to the tunnel barrier between layers. Parameters 
values can also be influenced by disorder on length scales 



shorter than those on which this coarse-grained contin- 
uum theory is applied; disorder on longer length scales 
would have to be treated explicitly as we mention in the 
discussion section. The upshot is that the numerical val- 
ues of the parameters in Eq.( [3]) are usually not accu- 
rately known and likely vary substantially from sample 
to sample. It is worth noting that n must vanish at finite 
temperatures when At — >■ since two-dimensional sys- 
tems cannot support spontaneous long-range phase or- 
der. Among all continuum model parameters the value 
of Ps, which is typically ^ 10~''eV, is likely the most 
reliably known. 

Since the which layer pseudospin and the true elec- 
tron spin have identical quantum mechanical descrip- 
tions, we can borrow from the ferromagnetic metal 
spintronics literatur o^^i^'^ and use the Landau-Lifshitz- 
Slonczewski(LLS) equation to describe how the semiclas- 
sical pseudospin dynamics is influenced by a transport 
current: 
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{j is the number current density for electrons.) The sec- 
ond term on the right- hand- side of Eq. ( S]) captures the 
transport current effect. Its role in these equations is 
similar to the role played by the leads in the pioneering 
analysis of coherent bilayer tunneling by Wen and Zee^i 
who argued that a term should be added to the global 
condensate equation of motion to account for the con- 
tribution of transport currents to the difference in popu- 
lation between top and bottom layers. Eq.(|4l) describes 
how a transport current alters the condensate equation of 
motion locally. Its justification for the bilayer quantum 
Hall case is discussed in more detail below. The essential 
validity of this expression has been verified by countless 
experiments. In the first term on the right hand side 



iJ^g = {2/hn){5E[m]/5m) 



(5) 



describes precessional pseudospin dynamics in an effec- 
tive magnetic field which is defined by the energy func- 
tional. The third term accounts for damping of the col- 
lective motion due to coupling to its environment, in the 
present case the Fermi sea of quasiparticles. The damp- 
ing term in Eq.(|31) has the standard isotropic form used 
in the magnetism literature. A microscopic theory^^ of 
damping in quantum Hall bilayers makes it clear that the 
damping is actually quite anisotropic. We return to this 
point below. 

Magnetic order in the metallic ferromagnets to which 
this equation is normally applied is extremely robust, 
justifying the assumption that the spin-density magni- 
tude is essentially unchanged even when the system is 
driven from equilibrium by a transport current. The 
only relevant degree of freedom is the spin-density direc- 
tion, whose dynamics is described by Eq.(|4]). We expect 
that the pseudospin transfer torque description of bilayer 
quantum Hall systems will be most reliable when the or- 
der is most firmly established, that is far away from the 
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phase boundary^' that separates ordered and disordered 
states. 

The transport-current (Slonczewski^^) term in Eq.(|4]) 
can be understood in several different ways. In the spin- 
tronics hterature this term is normahy motivated by an 
appeal to total spin conservation and referred to as the 
spin-transfer torque. The idea is that when a ferromag- 
net's spin-polarized quasiparticles carry a transport cur- 
rent through a spatial region with a non-coUinear mag- 
netization, they violate spin-conservation. The collective 
magnetization must therefore compensate by rotating at 
a constant rate which is proportional to the fermion drift 
velocity. The form we use for the pseudospin transfer 
torque assumes that each component of the pseudospin 
current is locally equal to the number current times the 
corresponding pseudospin direction cosine. This prop- 
erty does not hold locally in a microscopic theory, but is 
expected to be valid in the smooth pseudospin texture 
limit we address. (In metals an additional phenomeno- 
logical factor is required to account for the difference in 
drift velocity between majority and minority spin elec- 
trons.) 

Since pseudospin is not conserved in bilayer quantum 
Hall systems, as we can see from Eq.([2|) or Eq.([3l), this 
argument does not apply directly. If we appeal to a mean- 
field description of the pseudospin ferromagnet, however, 
we can obtain the same result by the following argument. 
Mean-field quasiparticles satisfy a single-particle Hamil- 
tonian for a particle in a magnetic field which experiences 
both a scalar potential, and a pseudospin-dependent po- 
tential that can be interpreted as a pseudospin effec- 
tive magnetic field. Following standard textbook deriva- 
tions it is possible^'^ to derive an expression for the time- 
dependence of the contribution of a single quasiparticle 
to any component of the pseudospin density. The ex- 
pression contains a pseudospin precession term and an 
additional term which is the divergence of the current of 
that component of pseudospin. Summing over all quasi- 
particle states we obtain a precession term that depends 
on the configuration of the order parameter, and hence 
on the pseudospin-field to which it gives rise, even in the 
absence of a current. We obtain the additional current- 
driven term on the right-hand side of Eq. ( |4]) when the 
pseudospin currents are non-zero and space dependent. 
The additional term in the equation of motion can also be 
viewed^^ as a consequence of an altered relationship be- 
tween pseudospin-dependent effective magnetic field and 
pseudospin polarization direction for quasiparticles that 
carry a current. Eq.(|3]) should in principle also include 
a current related damping term'^'^ which we ignore in the 
present paper. 

When quantum Hall bilayer pseudospin ferromagnets 
are tilted far from their easy plane, order tenda^ to be 
destroyed. For that reason we are often most interested in 
the limit in which rUz is much less than 1. It is therefore 
convenient to express the pseudospin direction in terms 
of the azimuthal angle ^, which is the inter-layer phase 
difference, and which is proportional to the layer po- 



larization. In terms of these variables the LLS equations 
take the form 
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In Eq.([n]) we have written j /n, which has units of veloc- 
ity, as the pseudospin velocity Vpg. 

These equations do not on their own provide a closed 
description of pseudospin dynamics in the presence of 
electrical bias potentials and need to be supplemented 
by a theory which specifies the spatial dependence of the 
pseudospin current. In general this quantity depends^! 
on the order parameter configuration as well as on the 
contact geometry and external current or voltage biases. 
The transport theory and the pseudospin dynamics the- 
ory are therefore not independent. In the present paper 
we study voltage biased Hall bars with source and drain 
leads at opposite ends, and with a variety of shapes and 
sizes. We assume, as a simplification, that the current 
distribution is defined by a local conductivity tensor with 
a large Hall angle. Given these simplifications, we are 
able to explicitly evaluate the maximum current at which 
time-independent order parameters are possible. Be- 
cause collective tunneling no longer contributes strongly 
to the dc interlayer current when the inter-layer phase is 
time-dependent, the interlayer conductance mechanism 
changes qualitatively when this maximum current is ex- 
ceeded. We therefore associate this current with the ex- 
perimental critical current. 



III. APPROXIMATE CRITICAL CURRENTS 

In this section we discuss approximate upper bounds 
on the critical current which are helpful in interpreting 
the numerical results described in the following section. 
We use a simplified version of the static limit of the ttIz 
LLS equation (Eq.( [5])) in which is assumed to be 
small: 

O^_|v2</>+i^sin0-ij.Vm,. (7) 

The three terms on the right-hand side can be identi- 
fied as contributions to the time-dependence of rriz (or 
equivalently of the exciton density) due respectively to 
the divergence of the exciton supercurrent, coherent con- 
densate tunneling, and the divergence of the z (layer 
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FIG. 1: (Color online) Separation of transport length scales 
in quantum Hall superfluid transport. This theory is intended 
to apply when the ordered state is well established and Hall 
angles are large because of the developing u = 1 quantum 
Hall effect. At large Hall angles current enters and leaves 
the samples at the hot spot corners, even when the source 
and drain contacts (gray) fully cover the ends of a Hall bar. 
When order is well established, the pseudospin orientation of 
the transport electrons achieves alignment with the conden- 
sate within a relatively small fraction of the sample area close 
to source and drain (solid yellow). In these areas pseudospin 
transfer torques convert transport currents into condensate 
counterflow supercurrents. When source and drain are con- 
nected to opposite layers a net supercurrent is injected into 
the interior region of the sample. In the remaining sample 
area (shaded blue) collective interlayer tunneling can act as a 
sink for the counterflow supercurrent. 



antisymmetric or counterflow) fermionic pseudospin cur- 
rent. The last contribution would be viewed as a spin- 
transfer torque in the analogous equations for an easy- 
plane anisotropy ferromagnetic metal. 

We start our quahtative discussion of critical currents 
by identifying some relevant length scales. First the 
length scale. 



A = 
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often referred to as the Josephson length because of the 
similarity between these equations and those which de- 
scribe Josephson junctions, emerges from balancing the 
first and second terms. In this paper we assume that 
the pseudospin magnetization direction departs from the 
x—y plane only over a small region close to the source and 
drain contacts whose spatial extent is small compared to 
the Josephson length. (This issue is addressed again in 
the discussion section.) If this is correct, we can separate 
length scales by identifying a region close to the contact 
which is small enough that we can ignore coherent tun- 
neling by setting Atn to 0, and large enough that we can 
assume that ruz is close to zero in the rest of the system. 
(See Fig. ([1]).) When Atn — )■ 0, Eq.([7]) simply expresses 
the conservation of the sum of the excitonic and quasi- 
particle counterflow currents. Integration of Eq.([7]) over 
the area close to the contact then simply describes con- 



version of quasiparticle counterflow current into conden- 
sate counterflow current. The total counterflow current 
emerging from the area near the source contact is half of 
the total current flowing into the system since = ±1 
in the contact and far away from the contact. 

In a tunnel geometry experiment the same counterflow 
current is generated near the source and drain contacts 
located at opposite dies of the sample, so the total coun- 
terflow supercurrent injected into the system is equal to 
the total number current flowing through the system. 

With this separation of length scales the quasiparticle 
(pseudospin transfer torque) term can be dropped in the 
remaining portion A of the total system area A. For 
static solutions of the LLS equations the condensate must 
satisfy an elliptic sine-Gordon equation inside A: 



-sm6 = 0. 



(9) 



When A cxD {Atn — > 0), this equation states that W^cj) 
is zero. It then follows from Green's theorem that no net 
counterflow supercurrent can flow into the area A. In the 
tunnel geometry that means that time- independent order 
parameter values cannot be maintained in the presence 
of a transport current unless Atn ^ 0. The maximum 
tunneling current that can flow through the system is 
particularly simple to determine in the small Atn, large 
A limit. When A is much larger than the system size the 
angle (j) cannot vary substantially over the system area. 
With this simplification the elliptic sine-Gordon equation 
can be integrated over the area A to obtain 



Ps 



V(t> ■ n = sin(0) 



(10) 



where P is the perimeter of A and n is proportional to 
the outward normal. The left hand side of Eq.f fTUt is the 
net supercurrent which flows out of the region A from its 
boundaries near the source and drain contacts, identifled 
above as the number current flowing through the system. 
Since the maximum value of | sin((/))| is 1, it follows that 
the maximum current consistent with a time-independent 
order parameter in this case is 



11 



eAAtn 
2h 



(11) 



Since the critical current in the small Atn limit is pro- 
portional to the area of the system we will refer to this 
quantity as the bulk critical current, as suggested by the 
notation used above. 

For larger Atn, A is no longer larger than the system 
size and it is not possible to maintain the maximum value 
of sin((/)) across the system. The LLS equation critical 
current in this regime depends on geometric details and 
we have not been able to obtain rigorous bounds. We can 
make a rough estimate by following an argument along 
the following lines. The elliptic sine- Gordon equation is 
very similar to the regular sine- Gordon equation in which 
second order derivatives with respect to time and position 
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appear with opposite signs. This 1+1 dimensional sine- 
Gordon equation appears as the Euler-Lagrange equa- 
tion of motion of a system with a Lagrangian with a 
kinetic-energy term proportional to ps(dt4>{x,t))'^ and a 
potential energy term proportional AfU cos{(j){x, t). Since 
total energy (integrated over position x) is conserved by 
this dynamics it follows that the variation of the typical 
value of ps(dt4>Y along the space-time boundary cannot 
be larger than ^ Afn. When this energy conservation 
condition is mapped from the regular sine-Gordon equa- 
tion to the (imaginary time) elliptic sine-Gordon equa- 
tion we can conclude that the typical value of V0-n along 
the boundary of A near the source contact cannot differ 
from the typical value of V(/) • n along the boundary near 
the drain contact by more than ~ Ajn. It follows that 
the current flowing through the system from source to 
drain should not be much larger than 



eWps 



(12) 



where W is the length of the contract region, or the width 
of a Hall bar assumed to be contacted at its edges. We 
will refer to as the edge critical current, since it is 
limited by the length of one edge. For stronger inter- 
layer coupling then critical current is expected to vary 
as (Ajn)^/^ once A is smaller than the Hall bar length. 
Finally we note that because of hot-spot effects in trans- 
port with large Hall angles the pseudospin transfer torque 
will act at the sample corners. Since the supercurrent is 
converted into coherent pseudospin precession over the 
length scale A, the effective size of the contact region will 
be ~ A when the Hall angle is large and A is smaller 
than W. We therefore estimate that the critical current 
is close to 



(13) 



independent of Atn, under these circumstances. We refer 
to this last critical current as the corner critical current. 
In the following section we compare numerical LLS crit- 
ical currents with these rough estimates. 



IV. MODEL CALCULATIONS 

The numerical calculations we describe below are sim- 
ilar to those carried out in micromagnetic descriptions of 
ferromagnetic metal spin-transfer torque physics, but are 
applied here to pseudospin transfer physics in condensed 
bilayers. We divide the system area into pixels within 
which the pseudospin magnetization is assumed to be 
constant. Except where noted we used square 101 x 10^ 
pixels where I is the magnetic length. We rewrite the 
spin-transfer torque term in Eq.( [5]) in the discretized 
form: 



where k labels the four neighboring sites, A is the pixel 
area, and a labels components of the magnetization ori- 
entation in the pixel of interest. /„ is the quasiparticle 
number current flow from neighbor site n into this pixel. 
In this article we estimate values of /„ by solving the 
resistor network model obtained by discretizing a contin- 
uum model with a local conductivity that includes a Hall 
component. 



{/} = [G]{V} 



(15) 



In 



[mk,a ~ {rrik ■ m)ma] (14) 



where [G] is the conductance matrix that describes the 
local conductivity including Hall components, {V} is a 
vector of local voltages in each pixel, and {/} is a vec- 
tor composed of currents that flow between pixels. Us- 
ing current conservation conditions the dimension of the 
matrix can be reduced to the pixel number N. Given 
the source and drain contact voltages, we can solve for 
the internal voltage and current distributions and for the 
variation of source and drain currents across the contacts. 
Note that quantum Hall physics conies into play through 
the Hall contributions to the conductivity matrix [G]. 
The local Hall conductivity was set to axy = e^/h, which 
is close to the appropriate value for v = I whether or 
not the Hall plateau is fully formed, and the longitudinal 
conductivity was set to 

c^xx = 0.05 exp(— m^/Vt^) (m • mL){m ■ rnn) /h (16) 

where tol.a is the magnetization orientation to the left 
and right of a boundary separating two-pixels. The Hall 
angle used in these calculations was therefore tan^^(20) 
over the largest part of the sample in which the pseu- 
dospin magnetization is close to collinear and planar. 
The results we report on are not sensitive to the Hall 
angle, provided that it is large. 

The current which flows into a pixel is assumed to have 
the same pseudospin polarization as the pixel from which 
it is incident. Currents entering or exiting from the con- 
tacts are assumed to have rUz = 1 for top layer contacts 
and rUz = — 1 for bottom layer contacts. In this way 
the pseudospin transfer torque and the Landau-Lifshitz 
precessional torque acting on each pixel's pseudospin can 
be evaluated. Note that the pseudospin torque depends 
not only on current paths, but also on the pseudospin 
magnetization configuration. 



A. Critical Current Identification 

As mentioned in the previous section, we identify the 
critical current as the circuit current value above which 
a time independent solution no longer exists. When the 
external current is small, the Gilbert damping term in 
the LLS equation relaxes the pseudospin magnetization 
into time-independent configurations. The behavior of 
the pseudospin as the current increases is partly analo- 
gous to the behavior of a damped pendulum driven by 
an increasingly strong torque. 
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FIG. 2: Magnetization change per time step 5m x A'^ vs. ap- 
plied voltage in a 500/ X 400/ system for 2-KfnAt = IO'^-Bq. 
This curve was obtained by changing the bias voltage by 
SV = 2.5 X 10~^ fiY at each time step. Eq = /d, the energy 
unit used in all our calculations has a typical experimental 
value ~ 7 meV. 5m develops large amplitude oscillation when 
the voltage exceeds the critical voltage Vc- 



We identify the critical current numerically by slowly 
increasing the driving voltage (by SV per time step) and 
monitoring the change in magnetization. To be more 
explicit, we examine 



StueeY^ \m,{V + 5V) - mi{V)\/N 



(17) 



The sum in Eg.f I17p is over all pixels. The critical cur- 
rent can be defined precisely as the current above which 
5m remains finite when 6V — >■ 0. In practice we choose 
a suitably small value of SV and examine the current or 
voltage dependence of Sm. We find that Sm increases 
dramatically and begins oscillating at a voltage we iden- 
tify as the critical voltage. (See Fig.( [5]).) An alter- 
native but more laborious method of obtaining critical 
currents is to sequentially examine the dynamics of the 
pseudospin magnetization at a series of fixed values of 
the applied voltage V. If the applied voltage is below 
its critical value, Sm will approach zero exponentially at 
large times. If the applied voltage is above its critical 
value Sm will not approach zero and usually exhibits an 
oscillatory time dependence. In our calculations we used 
the first approach to determine an approximate value of 
the critical voltage (and hence the critical current) and 
the second method to refine accuracy. 



B. 



Ic vs. At 



It is useful to start by briefly discussing the single- 
pixel limit of the calculation, which should apply ap- 
proximately to the case in which the Josephson length 
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FIG. 3: (Color online) Critical current vs. 2iT£^nAt in a 
400Z X 500 I system. The dash-dot (red), long-dash (blue) and 
short-dash (black) curves plot the values of the bulk (/g), edge 
(Ie) and corner (Iq) limited critical currents discussed in the 
text for this sample geometry. The square dots plot LLS equa- 
tion critical currents at a series of Atn values. All the calcula- 
tions in this paper were performed using pseudospin stiffness 
(exciton superfluid density) ps = 0.005i?o where Eq = /el 
is the energy unit. Currents are in units of Iq — 1% ~ e.ps/h. 
For the value of ps used in these calculations Jo — 8nA. 



is longer than the system size. The steady limit of the 
LLS equation for the z-component of the pseudospin is 



1 AtnApixci . 
= — sm ( 
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m,M) (18) 



where / is the charge current flowing through the system 
and mz.L and mz.R are the z component of pseudospin for 
the source and drain leads at the left and right ends of the 
sample. For the drag geometry (current entering and ex- 
iting from the same layer) m^^L — mz^R, there is no spin 
torque term in the single-pixel calculation, and the steady 
state equation can be satisfied by setting sin0 — 0. For 
the tunneling geometry mz,L = —mz^R, the maximum 
pseudospin torque that can be compensated by the tun- 
neling term is obtained by setting siiKp — > 1. We there- 
fore obtain 



pixel 



(19) 



This gives a linear dependence of Ic on the single-particle 
tunneling strength Aj. The numerical procedures de- 
scribed above accurately reproduce this simple result. 

As explained previously and discussed more fully later, 
we believe that the pseudospin transfer torque in most 
quantum Hall superfluid experiments acts in a small frac- 
tion of the system area. We therefore need to perform 
calculations with many pixels in order to represent a typ- 
ical measurement. Fig. ([3]) shows numerical critical cur- 
rent results for a fixed sample geometry as a function of 
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FIG. 4: Spatial distribution of the z component of the pseu- 
dospin transfer torque in a system with 500 x 400 1? area, 
2nl'^nAt = 10"^£;o, and I/T = 0.29. The pseudospin trans- 
fer torque in the model studied here acts mainly in the hot 
spot pixels. Eo and / are defined as in previous figures. 



At. In this figure Iq = = eps/h ^ 8nA is the unit 
of current and Eq = /el is the unit of energy. Typi- 
cal quantum Hall superfluid experiments are performed 
at a magnetic field of roughly 2.1 Tesla for which the 
magnetic length is 17.65 nm which gives Eq = 6.4 meV. 
In all the calculations reported on here we used a pixel 
area ^pixci = 10 x 10 and the the mean- field theory 
estimate^ {ps — 0.005ii^o) for the pseudospin stiffness. 
The calculations in Fig. ([31) are for 40 pixels in the width 
W direction and 50 pixels in the current L direction. For 
this system size , the Josephson length is approximately 
equal to L when 27rZ^Atn ^ 4 x 10~''ii'o- The numerical 
results illustrated in Fig. (jS]) show the crossovers from 

1 /2 

Aj-dependence, to Aj -dependence, to saturation as At 
increases that was anticipated in our qualitative discus- 
sion. At small Atn the critical current is reduced by a 
small fraction compared to the single pixel result in ac- 
cord with the Fig.(l). 

We now examine the ingredients which enter these 
numerical results in greater detail. According to the 
schematic Fig. (Il}, the pseudospin transfer torque acts 
only near the hot spots at which current enters and ex- 
its the sample. Fig. Q shows a typical numerical re- 
sults for the spatial distribution of the pseudospin trans- 
fer torque. In the present model, the area of the region 
in which transport current is converted into supercurrent 
depends on the pixel size and the pseudospin stiffness and 
anisotropy cocfhcients. In the tunneling geometry, coun- 
terflow supercurrent is generated near both source and 
drain hot spots and flows diagonally toward the center 
of the sample. In Fig. ([S|) we show a typical supercur- 
rent distribution for the tunneling geometry case. The 
corresponding distribution for an equivalent drag exper- 
iment (both contacts connected to the same layer) is 
illustrated in Fig. ([SJ. Since the elliptic sinc-Gordon 
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FIG. 5: Supercurrent distribution in system with area 500 x 
400 ^^ 27rZ^nAt = IQ'^Eo, and I/F = 0.29. This plot is for a 
tunneling geometry in which the source is a top layer contact 
and the drain is a bottom layer contact. Supercurrents are 
generated near both hot spots and fiow diagonally toward the 
sample center. Eo and I are defined as in previous figures. 
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FIG. 6: Supercurrent distribution in system with area 500 x 
400 27r/^nAt = 10"''£^o. This plot is for the drag geometry 
case in which the source and drain are both top layer contacts, 
but other parameters of the calculation are identical to those 
used in the preceding tunnel-geometry figure. Eo and I are 
defined as in previous figures. 



equation applies locally when the pseudospin transfer 
torque is negligible, it follows from Green's theorem and 
Eq. ( [7]) that in a steady state the total countcrflow cur- 
rent injected into the interior of the sample (which for the 
tunnel geometry equals the total charge current flowing 
through the system) must match the area integration of 
{\/2){Atn/h) svtvcj). For a drag geometry experiment the 
same integral should vanish. 

We now examine steady state condensate configura- 
tions at currents near the critical current for both small 



8 



cp ( rad ) 




FIG. 7: Pseudospin phase distribution in a system with area 
500x400 ^^ 27r/^nAt = X^'^Eq and IjT = 0.75. The Joseph- 
son length at this value of At is ~ 2500/. Note that is 
roughly constant through the system and that its value is 
close to 7r/2 because / is close to Jc. The units used here are 
the same as in previous figures. 




FIG. 8: Phase distribution of pseudospin in a system with 
area 500 x 400 t\ 2Trl'^nAt = 1.2 x lO'-'Eo and I/T = 0.97 . 
The Josephson length at this value of At is ~ 701. Note that 
a steady state is reached even though <j) varies considerably 
across the sample and has values larger than n/2. The units 
used here are the same as in previous figures. 



nAt and large nAj limits. In the former case our expecta- 
tion that sin should be nearly constant except near the 
hot spots is confirmed in Fig. ([7]). Collective tunneling 
acts sinks the injected counterflow supercurrent at a rate 
that is nearly constant across the sample area. The large 
nAt case is more complex. (f> varies over a large range and 
sin((/)) changes sign in different parts of the sample. Near 
the sample corners ip varies rapidly with position because 
of the large counterflow supercurrents. In the interior of 
the sample the phase changes less rapidly because the 



flovi' pattern spreads and because coherent tunneling is 
providing the required current sink. The critical current 
in the large nAt hmit, depends in a complex way on the 
geometry of the sample and on the spatial distribution of 
the pseudospin transfer torques. Nevertheless, the crit- 
ical current saturation we find in our numerical studies 
suggests that once the Josephson length is smaller com- 
pared to both the width and the length of a Hall bar 
sample with a large Hall angle, it is no longer relevant to 
the critical current value. 



C. 7c vs. System Geometry 

Finally, we briefly discuss critical current dependence 
on Hall bar dimensions at flxed nAt- In Fig. ([5]) we plot 
/■^ vs. Hall bar length in a series of model samples with a 
flxed single pixel width W — 101 and single-particle tun- 
neling amplitude 27rPnAt = IO^^Eq. For these param- 
eters the W is much smaller than the Josephson length. 
The critical current increases linearly with sample length 
and is therefore proportional to sample area until it sat- 
urates at L 1000 L The length at which the current 
saturates is a bit longer than the Josephson length and 
the value of the critical current is accordingly somewhat 
larger than the ID estimate^! /|,. The large L behavior 
is however consistent with the expectation that the crit- 
ical current should not increase with system length once 
L is substantially larger than the Josephson length A. In 
Fig. pO)) we plot Ic vs. system width W with the length 
flxed at 500 1 and the same At as above. The range of W 
covered is limited somewhat by numerical practicalities 
and goes from a width much smaller than the Josephson 
length to a width which is somewhat larger. Over this 
range deviations from the ID model in which the critical 
current is proportional to Hall bar width are small. 



V. DISCUSSION AND CONCLUSIONS 

We start our discussion by commenting briefly on some 
essential differences between the critical current of a 
Josephson junction and critical currents for coherent bi- 
layer tunneling. In a Josephson junction, current can 
flow without dissipation across a thin insulating layer 
that separates two superconductors. The difference in 
condensate phase across the junction (pj is normally zero 
in equilibrium but can be driven to a non-zero steady 
state value when biased by current flow Ij in the cir- 
cuit in which the junction is placed. For thick insulating 
layers the current is related to the phase difference by 

Ij = 75 sin(</.j), (20) 

where Ij is the junction's critical current. Ea.f [20|) should 
be compared with Eq. ( [7]) . The most obvious difference 
is the appearance of the lateral 2D coordinate in the co- 
herent bilayer case. In the Josephson junction case, the 
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FIG. 9: Critical current vs. system lengtli L for a narrow 
Hall bar with width W = 10 I. The single-particle tunneling 
amplitude 2'jTl^nAt = 10~^(i?o) corresponding to A 250Z 
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FIG. 10: Critical current v.s. system width W for length L — 
500 I). The single-particle tunneling amplitude 27rZ^nAt = 
lO^^'(Eo), corresponding to A ~ 250/ 

lateral dependence of usually plays no role unless an 
external magnetic field is present. In the coherent bilayer 
case, on the other hand, lateral translational invariance 
is always broken because the pseudospin transfer torques 
that ultimately drive the coherent tunneling current do 
not act uniformly across the system. 

A closer comparison is possible in the special case in 
which the pseudospin stiffness is large enough to inhibit 
lateral variation of (p. Integration of Eq.([7l) over the area 
then yields for the coherent bilayer 

Ibl = Ih sin(0) (21) 
where Ig is the bulk critical current discussed in the body 



of the paper and 

iBL^e |^[|v0+^jm,] -h. (22) 

is the injected counterflow current. The most essential 
difference between tunneling in coherent quantum Hall 
bilayers and Josephson junctions lies in the difference in 
physical content between the bias currents Ij and Ibl- 
In the case of a Josephson junction it is a bulk dissi- 
pationless supercurrent which fiows perpendicular to the 
plane of the junction. In the case of a coherent quan- 
tum Hall bilayer, it is the counterflow (z) component of 
the quasiparticle current. The counterflow component 
of the quasiparticle current is normally fully converted 
to condensate counterflow supercurrents by pseudospin 
transfer torques, as we have discussed at length. The 
voltage drop across a Josephson junction can be mea- 
sured and vanishes below the critical current. Because 
fermionic quasiparticle transport, transverse counterflow 
supercurrents, and collective interlayer tunneling are un- 
avoidably intertwined in the coherent bilayer tunneling 
case, there is no corresponding measurable voltage which 
vanishes below the critical current. Instead, the critical 
current is marked experimentally by an abrupt increase 
in measured resistances. 

Next we draw attention to the significance of the qual- 
itative difference in experiment between drag geometry 
and tunnel geometry transport measurements. In our 
theoretical picture is correct there should be no qualita- 
tive difference difference between voltage measurements 
in these two geometries at currents below the critical cur- 
rent, once coherence is well established. (Important dif- 
ferences between these measuring geometries can occur in 
the interesting regime close to the phase boundary where 
interlayer coherence may be established over a small frac- 
tion of the sample area.—) Pseudospin transfer torques in 
tunnel and drag geometry experiments with similar total 
circuit currents give rise to similar counterflow supercur- 
rents, as seen in Figs. (5) and (6). Excitonic conden- 
sates have a maximum local (critical) supercurrent den- 
sity that they can support, which in the case of quantum 
Hall bilayer exciton condensates has been estimated^ to 
be ^ lAm~^. For a typical sample width of 10~^m this 
corresponds to a critical current in the mA regime, or- 
ders of magnitude larger than the current levels typically 
employed in quantum Hall experiments. Even account- 
ing for possible corrections due to disorder, local critical 
currents are unlikely to be approached experimentally 
for either contact geometry. (If they were, we would ex- 
pect more similarity between the two experiments.) The 
critical currents which are important for typical tunnel- 
ing experiments are not local critical currents, but global 
critical currents which set an upper limit on the rate 
at which injected supercurrent can be sunk by collective 
tunneling. Since there is no net injected supercurrent 
in the drag geometry experiment, only local limits apply. 
In the tunnel-geometry experiment, the global limit must 
be satisfled, setting a critical current scale which can be 
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orders of magnitude smaller if At is small. Although the 
relationship anticipated here between tunneling geome- 
try and drag geometry transport has not been specifically 
tested experimentally, it appears to us that it is clearly 
consistent with published data. 

We now turn to a comparison of our theory with exper- 
iment, and in particular with the recent experiments of 
Tiemann et al^ who have systematically studied the de- 
pendence of the tunneling critical currents in their sam- 
ples on temperature and layer separation. The data of 
Tiemann et al. appear to be broadly consistent with 
that reported in earlier~ii^ but focus more on tunneling 
anomalies in the region of the phase diagram far from the 
coherent state phase boundary. Tiemann et al. find i) 
that the critical tunneling currents in their samples are 
proportional to system area, ii) that their typical value 
is ^ lOnAmm^'^, and that they saturate upon moving 
away from the coherent state phase boundary^! either 
by lowering temperature or by decreasing the ratio of 
layer separation to magnetic length. 

Tiemann et a/.'s finding that the critical current is pro- 
portional to area would be consistent with our analysis if 
the Josephson length was comparable to or longer than 
the ^ 10mm length of their long-thin Hall bars. When 
in the bulk-limited critical current regime, the critical 
current should be given by 

^ = = 3.1S[Tesla] x KAcW] x lO^Amm"^ (23) 

where At = 27rPAt n is the interlayer tunneling ampli- 
tude suitably renormalized by quantum and thermal fiuc- 
tuations and B is the magnetic field strength dX v = 1. 
Inserting the magnetic field strength, experimental criti- 
cal currents can be recovered by setting At 10^ ^"^eV. 
Therein lies the rub. Although values for At on this scale 
do imply Josephson lengths that are on the mm length 
scale and therefore consistent with bulk-limited critical 
currents in the samples studied by Tiemann et al., they 
are four or more orders of magnitude smaller than values 
expected on the basis of theoretical estimates. Suspicion 
that there is a fundamental discrepancy is established 
more convincingly, perhaps, by observing that these crit- 
ical currents also correspond to At values several orders 
of magnitude smaller than seemingly reliable experimen- 
tal estimates"^'^ based on the inter-layer tunneling conduc- 
tance in similar samples at zero-field. The key issue then 
in comparing critical current theories with experiment 
appears to be explaining why they are so small. 

We expect on general grounds that the experimental 
value of Atn = At/27rZ^ should be renormalized down- 
ward by both thermal and quantum fluctuations and 
by disorder. Indeed according to the familiar Mermin- 
Wagner theorem, the order parameter n must vanish for 
At — )• at finite temperatures. The importance of ther- 
mal fluctuations is strongly influenced by kgT/ps. If 
mean-field theory estimates? of ps can be trusted, the 
value of Ps in typical experimental samples should be 
^ 3 X lO^^eV and ksT/ps should therefore be less than 



0.1 at the lowest measurement temperatures. At these 
low temperatures thermal fluctuations alone appear to 
be insufficient'^^ to explain the discrepancy even though 
At/ksT is certainly small. The experimental finding 
that the critical current saturates at low temperatures 
supports this conclusion. Similarly, theoretical estimates 
that quantum fluctuation corrections to the order param- 
eter are not^^ large well away from the transition bound- 
ary are consistent with the experimental finding of satu- 
rating critical currents in this regime. It appears that an 
explanation for the small critical currents must be found 
in the disorder physics of quantum Hall superfluids. 

The analysis of tunneling critical currents presented in 
this paper can accommodate disorder implicitly through 
its influence on the parameters ps and Atn. Including 
disorder effects through renormalized coupling constants 
would be adequate if the characteristic length scales for 
disorder physics are smaller than characteristic length 
scales like A which are relevant for pseudospin transfer 
torques. In quantum Hall superfluids disorder may play a 
more essential role by nucleating charged merons^ (vor- 
tices). As we have explained, provided that the pseu- 
dospin transfer torque acts only close to the source and 
drain, the critical current is proportional to the integral 
of Atn over the area of the sample. This integral is re- 
duced to zero by a single undistorted meron located at 
the center of a large sample. It has in fact ofteni^"— i^"— 
been recognized that disorder induced vortices might play 
an essential role in many of the transport anomalies asso- 
ciated with bilayer coherence. It seems likely to us that 
this type of physics is very likely responsible for the small 
critical currents seen in experiment, but that existing the- 
ory is unable to account for this effect quantitatively. The 
current status of the subject calls for a detailed analysis 
of how they influence critical currents. On the experi- 
mental side, the importance of complex disorder-related 
pseudospin textures for critical current values could be 
reduced and the essential physics which limits critical 
currents revealed more clearly by studying samples with 
much larger bare values of At. 

Although we have attributed the substantial quanti- 
tative disagreement between pseudospin transfer torque 
theory and experiment to disorder-induced pseudospin 
textures and have suggested a strategy for achieving more 
quantitative tests of our theory, it is appropriate to step- 
back and reconsider other theoretical pictures that might 
be relevant to coherent bilayer tunneling experiments 
under some circumstances. For example, the version 
of the pseudospin-transfer torque theory applied here is 
based on the simplest possible assumption for the local 
pseudospin-polarization of the transport current, namely 
that the pseudospin current polarization simply follows 
the pseudospin density polarization. This assumption 
is certainly not generically correct, but its replacement 
requires more detailed knowledge of fermionic quasipar- 
ticle transport behavior. One approach is to assume that 
the quantum Hall effect establishes edge state transport 
and use experimental voltage probe data to infer— the 
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length scale on which the pseudospin polarization of in- 
jected electrons is relaxed, and therefore the length scale 
over which the pseudospin transfer torque acts. The 
advantage of this approach is that experimental data 
could be used to obtain the spatial distribution of pseu- 
dospin transfer torques. Examination of coherent bilayer 
transport data suggests^ that the torques sometimes act 
along most of the perimeter of the system, as assumed 
in a previousi^ attempt to understand coherent-bilayer 
tunneling data, and sometimes close to the source and 
drain the contacts as assumed here. The current model 
is most appropriate well away from the coherent state 
phase boundary as we have discussed. 

The present version of the pseudospin transport torque 
theory does not account for thermal or quantum fluctua- 
tions of the condensate, which are unimportant in metal- 
lic ferromagnets but might sometimes be important in 
coherent bilayers. Weak-coupling theoriesi^"— of bilayer 
tunneling do account for fluctuations, but treat At per- 



turbatively. These theories cannot account for the exis- 
tence of a critical currents and in practice assume that 
each layer has a separate well-defined local chemical po- 
tential. It is clear from published transport data that this 
assumption is not always valid, in particular that it is not 
valid in the portion of the phase diagram far away from 
the phase boundary on which the present paper focuses. 
Experimentally'^^ the critical current value decreases as 
the phase boundary of the coherent state is approached 
by increasing either temperature or the effective layer 
separation d/l. It would be interesting to attempt quan- 
titative tests of the predictions of weak-coupling theories 
in the portion of the phase diagram close to the phase 
boundary. 
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